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Abstract 

Schlieren flow visualizations of transverse oscilla- 
tions of jets submitted to an acoustic perturbation 
are analyzed in this paper. The aim is to estimate 
the shape and the position of the median line of 
the jet. Two methods for image processing are pro- 
posed, based on complementary approaches : inter 
image and intra image analysis. Synthesized im- 
ages of an oscillating jet are used to validate each 
method and compare their performances in the case 
of noisy pictures. Illustrations are then shown on 
real laminar and turbulent jets. The results ob- 
tained using both methods are very close, showing 
their reliability. Applications investigated in this 
paper are focused on the estimation of the convec- 
tion velocity of perturbations along the jet, and the 
influence of the Reynolds number and of the chan- 
nel geometry upstream the jet formation. 



1 Introduction 

In flute-like instruments, acoustic oscillation is gen- 
erated by the coupling of the unstable modes of 
an air jet and the acoustic modes of a pipe res- 
onator. For soft blowing conditions, the jet is lam- 
inar and the instability is usually described, fol- 
lowing Rayleigh's approach ([Ij), under assump- 
tions of infinite inviscid 2-D flow. Certainly this 
approach needs to be combined with ad-hoc as- 
sumptions to cope with the specific geometry of 



the problem (semi-infinite jet, acoustic transverse 
perturbation). Furthermore, because it is linear, 
this perturbation approach may be limited to small 
amplitudes of jet deflection. For many instruments 
like transverse flutes and popular music flutes, the 
sound power required calls for higher total jet flux. 
Reynolds numbers, based on jet thickness and ve- 
locity, up to 10000 have been estimated in trans- 
verse flutes. Under such conditions, the jet be- 
comes rapidly turbulent downstream from the flue 
exit. The development of instability induced by the 
acoustic perturbation on such jets has been stud- 
ied by Bechert ([2], [3]) and Fletcher (H, [5]). Both 
propose semi-empirical models to describe the jet 
oscillations. 

The intimate detail of the coupling between the 
acoustic perturbation and the jet instability remain 
an unsolved problem. It appears that for both lam- 
inar and turbulent jets, the lack of theoretical mod- 
els make experimental data useful to fit, calibrate 
and evaluate semi-empirical models ([6 ). 

Different researchers have studied experimentally 
jet oscillations under acoustical perturbations using 
hot-wire anemometry measurements ( [7 , [8 for ex- 
ample). On the other hand, flow visualization has 
been widely used as an inspiration for a physical 
analysis and description of the fluid mechanisms 
and aeroacoustics at work in flutes ([9j). 

In the present paper, we present tools for image 
analysis of flow visualizations to extract quantita- 
tive data on jet instability under acoustic pertur- 
bation. The development of instability on a jet is 



experimentally investigated using a simplified de- 
vice (see section [2|): the pipe resonator is removed 
from the flute and replaced by an acoustic field gen- 
erated by loudspeakers. Using the Schlieren tech- 
nique, flow visualizations have been carried out for 
a variety of Reynolds and Strouhal numbers. 

Two different image analysis techniques have 
been developed (see sections [3] and |4|) and they 
are discussed in detail in this paper. The ex- 
tracted data is analyzed in the framework of a har- 
monic perturbation that is convected downstream 
and grows in the shear layers of the jet (see section 

ED- 

Finally, the influence of the Reynolds number 
and the channel geometry on the jet behavior are 
presented as an application in section [6l 

2 Experiment description 

Experimental setup: a jet, created by blowing 
through a slit, is acoustically forced by two out-of- 
phase loudspeakers (see figure [T]). The speed of the 
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Figure 1: Experimental setup. Jet exit, speak- 
ers and two microphones (acoustic velocity estima- 
tion). 

jet is controled by varying the pressure in the cavity 
Pf just before the jet formation. The jet central 
speed at the channel exit Ub, can be estimated using 
Bernoulli's equation: 

= (1) 

V Po 

where po is the fluid density. 



Jet structure: two Reynolds numbers Re = 500 
and Re = 3000 were chosen {Re = U^h/v^ with 
h = 10~^m the channel height and u the kinematic 
viscosity). They correspond to two different oper- 
ating modes of the jet: Re = 500 gives a veloc- 
ity around 7.5 m/s which corresponds to a normal 
blowing condition for a recorder and assures a lam- 
inar behavior, while Re = 3000 corresponds to a 
velocity around 45 m/s producing a turbulent jet. 

The acoustic excitation frequency / is ranging 
from 70 to 1400 Hz (when Re = 500) and from 
1200 to 2900Hz (when Re = 3000) to cover a 
Strouhal {St) range from approximately 0.05 to 
1.15 {St — fh/Uh)' A fixed amplitude of the acous- 
tic velocity is chosen for each of the two experi- 
ments described in the present paper and a typical 
value is 0.5% of the jet velocit}ll|. The acoustic 
velocity is measured using a Microflown velocity 
sensor (low amplitudes) and a microphone doublet 
( high amplitudes) at a position near the unper- 
turbed jet path, but out of the flow visualization 
window. Preliminary experiment had been carried 
out to check the acoustic velocity homogeneity in 
the jet and sensor area. 

Schlieren jet images: Schlieren technique ([10], 
section 3.2.2 pl34) is used to observe the behav- 
ior of the oscillating jet: using an optical scheme, 
light phase shift crossing an inhomogeneous media 
is converted into light intensity (see figure [2|). We 
use CO2 jet traversing the air, whose mass density 
assures enough difference to allow Schlieren visu- 
alization (see figure [3]) while producing a behavior 
close to that found in real flute-like instruments. 
Sequential images of the jet are taken with a digi- 
tal camera synchronized with a stroboscopic lighte. 
Its frequency is set such that multiples of its fre- 
quency are slightly out of phase with the excita- 
tion frequency providing an aliased representation 
of the oscillation. Frequency of the camera was set 
to 14fps, exposure time to Ijis. Approximately 100 
images are taken covering two cycles of jet oscilla- 
tion. Images are captured in raw, black and white, 
bmp files with size 1280 x 448, and 8 intensity bits. 



-•^This is considerably lower that the 10% observed in real 
flutes but allows us to cover an interesting range of frequen- 
cies without over exciting the speakers. 
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Figure 2: Schlieren basic implementation: a subset 
S of the test area with different medium (here the 
CO2 jet) deviates Ught rays. A knife edge, placed in 
the focus of the second lens blocks part of the rays 
having crossed S, thus producing intensity varia- 
tions on the screen (see figure [3]). 
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Figure 4: Intensity level of a single column from a 
Schlieren image (left) and cross-correlation of two 
intensity graphs corresponding to the same column 
in two successive images of the sequence (right). 




Figure 3: Schlieren jet image (see figure [2] for the 
optical principle). 

3 Image processing 

The goal of the image processing algorithm is to 
detect transverse jet displacement for every image 
in the oscillation period. Two algorithms are pro- 
posed in this section. Their performances will be 
checked and compared in section [H and [5] respec- 
tively. 

3.1 Cross-correlation method 

A first algorithm based on cross-correlation be- 
tween successive images is proposed. 

Let us consider the intensity graph of one col- 
umn of the image as shown in figure [H (left picture). 
For that particular column, the shape of the graph 
does not change much as time goes on: the shape 
is only shifted, according to the transverse position 
T] in the y direction due to the movement of the jet. 
However, the particular shape varies with the col- 
umn considered since it is determined by the mass 
distribution of the jet as shown in figure [3l 

Cross-correlation is calculated between two in- 
tensity graphs corresponding to the same column 
of two successive images of the sequence (see right 
part of figure [4j). The position of the peak allows to 



measure the displacement of the jet between these 
two successive images. This is repeated for each col- 
umn of the image. Since changes on the jet position 
may be smaller than one pixel, a parabola is fitted 
to the three highest points to estimate the position 
of the peak. When considering all the columns of 
all the images in the sequence, the jet position can 
be reconstructed as shown in figure [5] {Re = 500 
and St = 0.2). 



Detaclad ampliti.da. luinere crorte. sf = 0.13305a. 
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Figure 5: Jet deviation, as estimated after image 
analysis through cross-correlation approach, in the 
case Re = 500 and St = 0.2 

Some defects of the experimental setup are 
clearly visible in left part of figure [H the two peaks 
corresponding to the brighter and darker part of 
the image should ideally be of the same height, but 
in practice a symetrical contrast is very difficult to 
obtain. Moreover, unhomogenous background in- 
tensity (inside each image) is represented as a non 
horizontal line on the borders of the left image in 
figure m However, these features of the images are 
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useful for the cross-correlation since they help to 
produce the well shaped peak in the curve observed 
in the right part of figure [H 

3.2 Morphological method 

The idea of this method is to apply morphologic 
functions on binarized images, in order to identify 
large spatial scales, corresponding to the jet main 
structure, while ignoring the smallest scales (vor- 
tices, background noise ...). This is done through 
the five following steps (also illustrated in figure [H): 

Contrast enhancement and homogenization 
between all the images of a sequence (fig- 
ure [GId): the contrast is first enhanced by his- 
togram equalization; the histogram of each image 
is streched to lie on the full 8-bit range, and to 
match a flat histogram. This also plays the role of 
contrast homogenisation between different images 
of a given sequence. 

Conversion from greyscale to black and 
white images (figure [6t) : a statistical method 
based on the Otsu principle estimates the thresh- 
old which best separates the histogram within two 
classes (pT]). The binarization threshold is there- 
fore different for each image. As a global method, 
the efficiency of the Otsu method is known to vary 
according to the type of image. A classical improve- 
ment is adopted, which consists in multiplying the 
Otsu threshold by an empirical constant in order 
to better enhance the jet structure for a given se- 
quence of image. Typical values for this constant 
lie beween 0.9 and 1.1. In the resulting BW im- 
ages, the jet is represented by a lower and a upper 
half jet, black and white respectively (figure[6t). In 
the following, both half jets are treated separately : 
the same operations are applied to the BW image 
and to its 1-complemented image. For clarity, only 
operations on the image containing the lower half 
jet are shown in figures [6]i, [6^ and[6]F. 

Extraction of largest spatial scales in the 
image ( figures [6]i and [6^): for each half-jet, 
the black and white picture includes several black 
structures, among which the half jet. This image 
is usually noisy, with black dots or small-sized sets 
of black pixels on the background. This can result 



either from a non-optimal threshold, or from the 
Schlieren experiment itself, which does not produce 
perfectly uniform backgrounds (mainly because the 
optical lens are never dispersion free). Therefore, a 
morphological closure ([12j) is performed to clean 
all small structures (figure [6li). Then a morpholog- 
ical opening is applied in order to highlight large 
scales in the image (figure [6^). The structuring el- 
ement is a disk, the size of which (optimized man- 
ually), is the same for all the images of a given 
sequence. 

Identification of the half-jet among the re- 
maining large scales (figure [6f): The image 
is decomposed into different regions, defined as 
groups of contiguous pixels. A unique region, cor- 
responding to the estimated half-jet, must be re- 
tained. However, since the jet is oscillating, some- 
times with large amplitudes, its position may vary 
significantly between different images of the same 
sequence. Therefore the discriminating criterion re- 
tained, involves both the area and the localization 
of each region and uses the (known) main flow di- 
rection from left to right). 

Jet edges and median line estimation : jet 

edges detection in figure [6t is trivial and is per- 
formed with Sobel algorithm which uses gradient 
information ([l3]). Edge detection is interesting in 
order to analyse in further studies the spreading of 
the jet as a function of the distance from the exit. 
Moreover, it is an intermediate step to estimate the 
median line, which we are mainly interested in. Fi- 
nally, an image containing edges of both half jets is 
scanned column after column and the median line 
is constructed point after point, including empirical 
criteria to select, among several possible paths, the 
one to follow. This guarantees that the resulting 
median line is single- valued, but possibly discon- 
tinuous in the case of vortices. 

Structure length: The largest the distance from 
the exit, the most complex the structure of the jet, 
and the least reliable the estimation of the median 
line. An index of confidence is then used: the es- 
timation of the median line stops at the distance 
(i.e. the column) corresponding to the length of 
the smallest half-jet. Indeed beyond that limit, 
median line estimation would rely on information 
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(6c) Automatic binarization (6d) Morpho. closure 




(6e) Morpho. opening (6f) Jet identification 




Figure 6: Morphological method: illustration of the 
different steps (from left to right and top to bot- 
tom). 

given by a single half-jet (as can be seen in fig- 
ure [7j). This distance depends on the experimental 
conditions (Reynolds and Strouhal numbers) but 
obviously varies from an image to another within a 
sequence, depending both on the phase of the exci- 
tation (whether the jet is horizontal or defiected), 
and on the image itself (the whole process of jet 
identification being more or less efficient from an 
image to another). 

The result of this method on one image is pre- 
sented in figure [71 Bringing together the me- 
dian line found for each image allows to construct 
3D plots similar to those obtained through cross- 
correlation image analysis (figure [5]) . 

4 Evaluation of both methods 



Both methods are evaluated in this section: they 
are first validated through the use of computer 
made images (section 14. ip and their robustness 
against background noise is compared in section 



Figure 7: Original Schlieren image superimposed 
with edge detected (black line) and the median line 
deduced (white line). 

4.1 Validation with digitally synthe- 
sized jet movies 

To validate both methods, the image analysis pro- 
cedures described in section [3] are applied to digi- 
tally synthesized images of a jet (see figure [8j). 

The images have been synthesized to imitate 
Schlieren images of real laminar jets like the ones 
captured in our experiments: same lines and 
columns numbers, 8-bit intensity images, upper 
and lower shear layers of the jet respectively black 
and white, grey homogenous background. 

The synthesized jet has an exponentially ampli- 
fied sinusoidal motion with controlable parameters. 
It is checked for different parameters of the jet (jet 
width, amplification coefficent, frequency of oscil- 
lation) that both methods are sucessful (see for ex- 
ample a graphical confirmation in the first row of 
figure[8]where the median line is superimposed with 
the synthesized image). 

4.2 Sensitivity to background noise 

Since a real image is always noisy, it is worth test- 
ing the infiuence of noise on the image analysis re- 
sults. To do this, white noise has progressively been 
added to the original synthesized image. This gives 
some insight on how algorithms behave when back- 
ground noise is increasing and contrast is decreas- 
ing. Results are presented in figure [8] for different 
signal to noise ratio. 



(8e) -3dB noise added (morphological analysis) 



(8f) -3dB noise added (crosscorelation analysis) 




(8g) -OdB noise added (morphological analysis) 



(8h) -OdB noise added (crosscorelation analysis) 




(8i) +3dB noise added (morphological analysis) 
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(8j) +3dB noise added (crosscorelation analysis) 



The morphological approach (left column) ap- 
pears to estimate successfully the median line on 
the whole image length until —3dB of noise. When 
the noise ratio is increased to OdB^ the median line 
is still estimated but not until the right border of 
the image. In the noisiest case {-\-3dB) the range 
on which the median line would have been esti- 
mated would be even shorter. To show how to pre- 
vent this, we have slightly modified the parameters 
of the analysis: the empirical constant multiplying 
the Otsu binarisation threshold is increased up to 
1.5 for the upper half jet and decreased down to 
0.5 for the lower half jet. Moreover, the radius of 
the structuring element for the closure has been di- 
vided by two. This allows to track the median line 
on the whole length by manual adjustment of the 
analysis parameters, but the estimation is somehow 
more crude. 

The crosscorrelation method gives reliable results 
for moderate levels of noise but shows more sensi- 
tivity than the morphological method for the nois- 
iest images. However, real Schlieren images ob- 
tained in our experiment are less noisy than the 
extreme case tested here. 

5 Results comparison 
5.1 Data analysis 

Every column of the images oscillates in time at 
the same frequency of the excitation. It is therefore 
possible to fit a sinusoid of that frequency to each 
column, and obtain the amplitude and phase 

Zy(/) of the fitted jet transverse displacement: 

where fs is the sampling frequency, N is the num- 
ber of images in the sequence, is the jet deflec- 
tion for that column at image number k. These 
values Y{f) are used to recreate a fitted version of 
the jet position as shown in figure [9l 

In a linear description of the jet ([14J, [15j), 
the perturbation is amplified while being convected 
downstream. The perturbation travels at a veloc- 
ity about one half of the jet speed. Inspired from 
this theory, the experimental data can be analyzed 
assuming a jet transverse displacement rj following 



Rtted amplitude + ^unsealed: awusNc speed, lumiere droits, str = 0. = 23S. 250000 




tim6[£6L:] " " Cistansefrcmthe flje exit [mm] 

Figure 9: Fitting of a sinusoidal model (cf. Eq. 
([2])) on the data plotted in figure [5l The wave su- 
perimposed at the flue exit corresponds to the per- 
turbating acoustic velocity. It is plotted for phase 
reference. 

my- 

r]{x, t)=^ (r/oe^^e''^^'"^/"^^) (3) 

where r]o represents a complex initial amplitude of 
the oscillation, 7 is the spatial growing rate, Cp is 
the convection velocity and x is the distance from 
the flue exit. 

The growing rate 7 can be estimated by fitting an 
exponential to the detected amplitude curve. The 
fitting is performed on a range of distances from 
the flue exit which is chosen as follows: the start- 
ing point corresponds to the minimal resolution ob- 
tained with the two methods (that is half a pixel for 
the morphological approach) whereas the last point 
considered is such that the oscillation amplitude 
reaches 80% of the maximum oscillation amplitude 
estimated by the morphological method. This up- 
per limit was found empirically to correspond to the 
limit above which an exponentially growing model 
is no more relevant. 

The convection velocity Cp of the perturbation 
can be estimated as the slope of the linear fitting 
to the phase (on the same fitting range as for the 
estimation of the growing rate 7) as shown in fig- 
ure [To] for the cross correlation method and the 
morphological method, in the case of a laminar 
jet (Re=500, left column), and for a turbulent jet 
(Re=3000, right column). 
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Amplitude of the oscillation , str = 0.199958,f = 238.280000 
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Amplitude of the oscillation , str = 0.210473,f = 1568.300000 
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Figure 10: Amplitude (top) and phase (bottom) of the jet transverse displacement using the cross 
correlation method (pink line) and the morphological method (green line). In the left column, the jet is 
laminar (Re=500, excitation frequency 238 Hz, i.e. St=0.195). In the right column, the jet is turbulent 
(Re=3000, excitation frequency is 1568 Hz i.e. St=0.21). 



5.2 Methods comparison 

5.2.1 Laminar jet 

A laminar jet (Re=500) is first considered. We dis- 
play in the left column of figure [10] amplitude and 
phase of the sinusoidal fit of the jet displacement. 
The two methods show very similar results: 

• the oscillation amplitude follows roughly an 
exponential up to a distance oixjh^W where 
h is the thickness of the slit. 

• the phase shift of the oscillation decreases 
roughly linearly downstream from fine exit, 
corresponding to the delay induced by the con- 
vection of the perturbation on the jet at a 
somehow constant velocity. 

5.2.2 Turbulent jet 

The case of a turbulent jet (Re=3000) is now inves- 
tigated, for a similar Strouhal number [St = 0.21). 
Here again, the two methods give similar results 
(see figure [TOl right column). It should be pointed 
out that, due to the lower relative acoustic exci- 
tation amplitude in the turbulent case {vac/Uj = 
2.710"^, compared to Vac/Uj = 510"^ in the 
laminar case), the transverse jet displacement is 
much smaller on the right column of figure [TOl than 
on the left column {rjmax/h ~ 0.3 compared to 



Vmax/h ~ 0.8). Analyses in the right column of 
figure [To] therefore correspond to a very severe case 
of a turbulent jet with very small transverse dis- 
placement. Regarding the phase shifts, a down- 
ward bending is observed, as excepted in the case 
of a rapidly slowing turbulent jet. 



5.2.3 Singular flow behavior 

For some conditions, the fiow instability moves 
away from the simple exponential model described 
in the previous section. For laminar condition, this 
is the case for the roll-up of vortices as well as for 
the development of varicose modes on the jet. 



Roll- up of the jet: In the case of strong jet am- 
plification, the accumulation of the vorticity in the 
jet shear layers at the infiexion points induces the 
formation of vortices: as a result, the jet breaks 
into an alternated vortex street ([l7]). Depending 
on the Reynolds number, this vortex street may 
rapidly be decomposed by turbulence. Even if the 
cross correlation appears to track the jet transverse 
displacement further downstream than the mor- 
phological method, figure [TT] (left column) shows 
that both analysis methods used are affected by 
this vortex formation. 
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Developpement of varicose modes: Depend- 
ing on the symetrical properties of perturbations 
on both shear layers of the jet, the jet oscihation 
can be sinuous or varicose. Sinuous oscillations are 
characterised by a transverse displacement of the 
jet while varicose oscillations correspond to an os- 
cillation of the jet thickness ([15j). Even if sinuous 
modes appear to be dominant ( J8]) in flute-like in- 
struments, varicose contributions can sometime be 
observed related to fluctuations in the jet velocity. 
In the case of the presence of varicose oscillation 
on the jet (see figure \TT\ right column), the mor- 
phological method appears to track the jet median 
line further downstream than the cross-correlation 
method. 

A partial conclusion of these examples above is 
that, in both cases, the bias between results of the 
two methods is an indication for singular flow be- 
havior. 

6 Selected results on the jet 
behavior 

Two applications of the proposed methods are pre- 
sented within this section. 

6.1 Influence of the Reynolds num- 
ber on the convection velocity 

The results obtained with cross-correlation and 
morphological methods can be compared in a more 
general way by considering the evolution of the es- 
timated parameters as function of the frequency in 
the case of a laminar (figure [T2|) and a turbulent 
(figure [T3|) jet. The mean convection velocity of 
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Figure 12: Convection velocity of the perturbations 
on the jet as estimated through cross-correlation 
and morphological methods, laminar jet (Re=500). 
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Figure 13: Convection velocity of the pertur- 
bations on the jet as estimated through cross- 
correlation and morphological methods, turbulent 
jet (Re=3000). 



the perturbation on the jet (defined as the mean 
slope of the phase shift) appears to be nearly in- 
dependent of the frequency of oscillation. This is 
true for a laminar jet as well as for a turbulent jet. 
In the former case, Cp/Ub 0.4 and in the latter 
Cp/Ub — 0.3. These results have to be compared 
with the Rayleigh ideal linear jet for which Cp/Ub 
varies between 0.3 and 0.5 depending on the jet 
velocity profile. The values found are also in agree- 
ment with values suggested by other authors ([15]. 

L7J)- 

A closer look at the evolution of the phase with 
the distance x from flue exit indicates a higher con- 
vection velocity at short distances from flue exit (a 
downward bending is observed). This can be in- 
terpreted both in terms of the slowing down of the 
jet velocity and of the smoothing of the jet veloc- 
ity profile, both induced by the viscosity. A fit of 
the phase with a power law reveals that the speed 
of the jet decreases as a function of x~^-^^ in the 
laminar case (mean value obtained from all the ex- 
citation frequencies available) and as a function of 
^-0.53 ^Yie turbulent case (mean value on all ex- 
citations frequencies also). From [19j (pl33 and 
p323) the jet centerline velocity in the laminar case 
is expected to decrease as while the decrease 

is much faster in the case of a turbulent jet since 
the jet centerline velocity follows The fitted 

values appear to be very close to the theoretical val- 
ues, even if it should be noted that the dispersion 
around the mean value is important, due to exper- 
imental problems for some frequencies. However, 
this approach appears promising and could be ex- 
tended to tackle the detection of the transition from 
laminar to turbulence in a jet. 
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Amplitude of the oscillation , str = 0.367354,f = 448.280000 
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Figure 11: Two illustrations of singular flow behaviours (one image of the sequence (top) and analysis 
results). Left column: roll- up of the jet (Re=500, St=0.36). Right column: development of varicose 
modes on the jet (Re=500, St=0.68). 

6.2 Influence of the channel geome- 



Propagation velocity 



try on the convection velocity 



Another application of the methods presented can 
be illustrated by the comparison of different flue 
channel geometries, as studied by Segoufin [20j. 
Convection velocities of perturbations on jets is- 
suing from long (27mm) and short (1mm) channels 
are compared. Both channels show the same thick- 
ness h = 1mm. In the case of the short channel 
with squared exit, the convection velocity is clearly 
higher than in the other cases, confirming the ob- 
servations done by Segoufin [20j in the context of 
edge-tones. This can be interpreted, as suggested 
by Segoufin, in terms of the boundary layers' thick- 
ness of the jet, which becomes thinner in the case 
of short channel. The present experiment was car- 
ried using the channels also used by Segoufin, and 
our results strongly correlate with Segoufin's results 
and interpretation. 



short 

45 degree 

rounded 

squared 



^1 A 



+0 +0 
( A A 



Figure 14: Dimensionless velocity, flue exit short 
(dashed), squared (dashdot), rounded (solid) and 
45 degrees (dotted). 



10 



7 Conclusion 

Compared to hot-wire anemometry, flow visualiza- 
tion allows a non intrusive measurement as well as 
a clear identification of the jet as composed of fluid 
particles issuing from the flue channed (streaklines, 
[21j). On the other hand, flow visualizations does 
not give direct information on the velocity field 
(streamlines) and requires the use of a gas of dif- 
ferent light refraction than that of the surrounding 
air. 

The images obtained by the Schlieren method 
are extremely sensible to variations on the setup. 
The position of the blade, the aperture of the di- 
aphragm, the alignment of the optics, the satura- 
tion curve of the acquisition software, the vertical 
and horizontal position of the jet, etc. are some of 
the parameters that needs to be adjusted for each 
experiment. Small differences on these parameters 
produce image variations. Therefore, even under a 
attentive care, images obtained from any two ex- 
periments are different. 

Both methods allow to estimate fairly well the 
amplitude as well as phase evolution of the per- 
turbation. The main difference is that the cross- 
correlation method performs an inter-image anal- 
ysis while the morphological method is an intra- 
image technique. Moreover, to achieve its best re- 
sults the morphological method needs to be man- 
ually calibrated for each sequence of images while 
the cross-correlation method can be completely au- 
tomate. The resolution of the detection is restricted 
to one pixel in the morphological method while 
continuous for the cross-correlation. Turbulent jet 
and noisy images seem to be better handled by the 
morphological method while the cross-correlation is 
preferred for laminar jets. 

The slowing down of the jet can be quite accu- 
rately evaluated, and appears to follow the theoret- 
ical values expected for both laminar and turbulent 
jets. 

The two methods we developped give very close 
results in terms of jet amplification and convec- 
tion velocities for a wide range of control param- 
eters. Together with the interpretation of results 
presented by Segoufin using the same channels, this 
indicates that both methods can be reliable and 
useful to the study of jet oscillation in edge-tones 
and flute-like instruments. 
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